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During  the  past  year  particularly  striking  progress  has  been  made 
toward  the  dual  objectives  of  this  program:  1)  to  measure  absolute  rate 
constant  temperature  dependences  for  radical-molecule  reactions  of  specific 
importance  to  modelling  of  combustion  and  plume  processes  in  aircraft  and 
rockets;  and  2)  to  contribute  to  the  general  understanding  of  the  temperature 
dependence  of  bimolecular  reactions,  particularly  as  related  to  observed 
curvature  in  Arrhenius  plots.  Improved  sensitivity  in  the  flash  photolysis- 
resonance  fluorescence  technique  was  achieved  and  high-accuracy  rate  con- 
stant measurements  were  made  in  the  heretofore  only  marginally  accessible 
500-1000K  temperature  regime.  Kinetic  data  for  the  important  reactions 
OH  + H2,  CH*,  CO,  and  C2H6  was  obtained  and  promising  future  experiments 
were  thereby  defined. 
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Aircraft  and  rocket  exhaust  plumes  provide  an  energetic,  highly  non- 
uniform,  local  environment  within  which  microscopic  physical  and  chemical 
processes  act  to  dictate  the  nature  of  bulk  observations  of  the  medium. 

Analyses  of  fuel  performance  and  tracking  procedures  rely  on  modelling  of 
those  contributing  processes  which,  in  combination,  govern  the  efficiency 
of  the  propulsion  system  or  surveillance  method  under  examination.  Chem- 
ically, reactions  within  the  combustor  and  afterburning  within  the  plume 
place  heavy  demands  on  the  estimating  prowess  of  kinetic  modellers.  In 
most  cases,  it  is  simply  not  possible  to  calculate  which  among  the  myriad 
of  possible  processes  are  most  Important  (and  thus  which  among  the  possible 
molecular  species  are  of  greatest  concern)— an  extensive  data  base,  promot- 
ing realistic  modelling  parameter  limitations,  is  needed.  From  a kinetic 
standpoint,  optimized  acquisition  of  this  data  base  Implies  the  measure- 
ment of  oxidant/fuel /reaction  Intermediate  reaction  rate  constants  needed 
directly  for  model  input  and/or  required  as  bases  for  the  generation  of 
'global'  kinetic  equations  applicable  to  entire  families  of  chemical 
reactions.  Improvement  of  the  understanding  of  those  factors  which  govern 
the  temperature  dependence  of  combustor/plume  relevant  blmolecular  reac- 
tions is  the  primary  goal  of  this  research  program. 

PROGRESS  SUMMARY 

Historically  speaking,  our  modern  view  of  the  temperature  dependence 
of  the  rate  constant  for  blmolecular  gas-phase  reactions  originated  with 
Arrhenius'  model  In  1889.*  Arrhenius'  postulated  equation,  k(T)  = Aexp-(E  /RT) , 

a 

endured  as  law  Into  the  1960s.  The  development  of  improved  experimental 

__i  1 
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techniques  applicable  over  extended  temperature  ranges  of  measurement  then 
spawned  a large  number  of  intercomparable  studies  and  led  to  the  initial 
clear  case  of  failure  of  Arrhenius'  equation.  Even  up  through  the  past 
few  years,  however,  such  failures  were  typically  characterized  merely  by 
noting,  for  a particular  reaction,  substantial  disagreement  between  the 
log  k (T ) vs.  1/T  slopes  obtained  as  the  high- temperature  (T  > 1000K)  and 
low- temperature  (T  < 500K)  straight-line  extrapolations.  It  has  not  been, 
in  the  past  decade,  that  curvature  in  the  log  k(T)  vs.  1/T  (so-called 
Arrhenius)  plot  is  unacceptable;  theoretically,  curvature  has  been  con- 
sidered likely  for  some  time.  The  lack  of  sufficient  congruence  of  data 
in  either  (high  or  low)  temperature  regime  and  the  nearly  total  absence 
of  reliable  rate  constant  data  between  500K  and  1000K  has  delayed  until 
quite  recently  and  limited  to  just  a very  few  reactions  even  the  empirical 
smooth  joining  of  high  and  low  temperature  measurements. 

It  is  precisely  that  void  In  elementary  reaction  rate  constant  measure- 
ments from  500K  - 1000K  that  has  Inspired  the  development  and  growth  of 
this  research  program.  During  the  past  year,  accurate  extended  temperature 
range  rate  constant  measurements  have  been  made  for  the  following  reactions: 


OH 

+ H2  - 

ki  y 

HjO 

+ H 

(1) 

OH 

+ CH„  - 

HjO 

♦ CH, 

(2) 

OH 

+ co  - 

H 

+ COj 

(3) 

OH 

+ CjH*  • 

HjO 

♦ c2h$ 

(4) 

Distinct  curvature  in  the  Arrhenius  plot,  varying  In  degree,  has  been 
observed  for  each  Reaction  (1)  - (4). 
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Our  results  for  Reaction  (1)  are  listed  in  Table  I and  are  plotted 
individually  and  in  combination  with  selected  other  results  In  Figures  1 
and  2,  respectively.  Our  measurements  demonstrate  that  the  degree  of 
curvature  in  the  Arrhenius  plot  for  Reaction  (1)  exceeds  that  predicted 
by  the  expression  of  Cohen  and  Westberg.  Cohen  and  Westberg's  empirical 
fit  was  generated  primarily  by  interpolative  joining  of  immense  numbers  of 
measurements  (not  shown  in  Figure  2)  in  the  low  and  high  temperature 
regimes.  Thus,  even  though  the  kinetics  of  Reaction  (1)  has  received 
enormous  attention  over  the  years,  the  previous  scarcity  of  data  from  500- 
1000K  has  limited  the  reliability  of  the  best  obtainable  empirical  fits 
for  MT). 

Absolute  reaction  rate  constant  measurements  obtained  for  Reaction  (2) 

are  listed  in  Table  II  and  plotted  In  Figures  3 and  4.  Many  Important 

general  features  are  evident  from  Figure  4.  First,  it  can  be  seen  that 

straight-line  extrapolations  of  either  the  low-temperature  based  (NASA, 

troposphere  and  stratosphere)  recommendation  or  the  high-temperature 

based  (Peeters  and  Mahnen  best  fit)  recommendation*  produce  order  of 

magnitude  errors  in  k2(T)  In  the  opposing  temperature  regimes.  Clearly, 

even  for  approximate  modelling  efforts,  curvature  in  the  Arrhenius  plot 

cannot  be  ignored.  Second,  it  is  interesting  to  contrast  the  data  of 

c 

Zellner  and  Stelnert  with  Zellner's  subsequently  obtained  forced  empirical 
fit®  to  the  existing  high-  and  low-temperature  data  (solid  line  In  Figure  4). 
Using  the  direct  experimental  technique  of  flash  photolysis-resonance 
absorption,  Zellner  and  Stelnert  measured  k2(T)  up  to  temperatures  of  892K. 

At  lower  temperatures,  this  data  agrees  well  with  the  empirical  fit,  but, 
at  higher  temperatures,  their  measured  k2(T)  values  exceed  the  best  fit 
values  by  a factor  increasing  with  temperature  and  already  more  than  two 
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Table  I 


oh  + h2-+h2c>h 


( 0 K ) 

TEMPERATURE 


( C M ^ MOLEC*^  S E C " ^ ) 

K B I M 0 L E C U L A R* 


2 9 8 

( 6 . 

A 9 9 

( 1 . 

5 7 6 

( 2 . 

6 A 8 

( 3 . 

7 3 9 

( 7 . 

8 3 8 

( 1 , 

9 0 A 

( 1 . 

9 9 2 

( 1 . 

0 8 

+ 

.37) 

X 

1 0 " 

1 5 

1 5 

+ 

.08) 

X 

1 0 ■ 

1 3 

3 A 

+ 

.18) 

X 

1 0 ’ 

1 3 

8 A 

+ 

.18) 

X 

1 o • 

1 3 

6 A 

+ 

.56) 

X 

1 0 ' 

1 3 

3 0 

+ 

.10) 

X 

1 0 " 

1 2 

8 6 

+ 

.35) 

X 

1 0 ' 

1 2 

9 9 

+ 

.19) 

X 

1 o • 

1 2 

* Stated  error  bounds  represent  the  2°  values- 
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$ Present  Work.,  2°  error  bds 
a Westenberg  and  deHaas 
Cohen  and  Westberg 


Figure  2 


IA3Lt  11 


oh+ch4-*  h2o+ch3 


( 0 K ) 

TEMPERATURE 


( C M ^ HOLEC*  ^ S E C * ^ ) 
KBIMOLECULAR* 


298 

(7.50 

+ 

. 6 ) 

X 

10‘15 

398 

(9.73 

+ 

.95) 

X 

10-1“ 

9 9 8 

(8.07 

+ 

1.1) 

X 

10-1“ 

511 

(1.95 

+ 

. 12  ) 

X 

1 0 " 1 ? 

529 

(1.67 

+ 

.06) 

X 

1 0- 1 3 

600 

(3.19 

+ 

.90) 

X 

1 0 _ 1 3 

619 

(2.75 

+ 

. 99  ) 

X 

10-13 

696 

(5.78 

+ 

. 58  ) 

X 

1 0 * 1 3 

772 

(8.37 

+ 

1.5) 

X 

10-13 

915 

(1.50 

+ 

.15) 

X 

1 0 " 1 2 

1020 

(2.0 

+ 

. 2 ) 

X 

lO-l2 

Stated  error  bounds  represent  the  2°  values. 


at  900  K;  indeed,  Zellner  was  forced  to  attach  minimal  weighting  to  his 
own  measurements  to  attain  agreement  with  high-temperature  measurements 
such  as  those  of  Peeters  and  Mahnen.  In  contrast,  our  k2(T)  values  reported 
above  lie  dead-on  Zellner's  empirical  expression,  verifying  both  the 
adequacy  of  the  functional  form,  k ( T ) = ATnexp-(E/RT),  and  the  basic  sound- 

a 

ness  of  previous  measurements  in  both  temperature  regimes.  Our  flash  photol- 
ysis-resonance fluorescence  technique,  then,  appears  more  appropriate  to 
higher  temperature  investigations  than  the  flash  photolysis-resonance 
absorption  method  of  Zellner  and  Steinert.  Although  both  experiments 
measure  the  disappearance  of  OH  radicals  in  real  time,  the  absorption 
technique  monitors  the  di fference  between  two  large  numbers  [Ot_0  ~ 1^)/ 

It  ],  and  thus  requires  larger  OH  concentrations  to  match  the  statistics 
obtainable  from  the  fluorescence  technique.  Such  larger  radical  concentra- 
tions enhance  the  probability  of  OH  disappearance  due  to  processes  other  than 
Reaction  (2). 

Experimental  measurements  of  k3 (T ) for  Reaction  (3)  are  listed  in  Table 
III  and  plotted  in  Figure  5.  The  OH  + CO  H + C0£  reaction  is  of  extreme 
importance.  In  combustion  processes.  Reaction  (3)  represents  the  principal 
mode  of  conversion  of  CO  to  CO^;  since  this  reaction  is  very  exothermic,  its 
kinetics  considerably  impact  the  time  profile  for  combustive  energy  release, 
and  thereby,  factors  such  as  ignition,  acceleration  and  propagation  of  flames. 
Also  in  the  atmosphere,  Reaction  (3)  is  the  principal  CO  -+■  C0,>  conversion 
pathway;  its  kinetics  and  possible  pressure  dependence  thus  impact  numerous 
Important  atmospheric  chemicals  cycles  (e.g.,  NO  •*  NO,,).  The  pictorial 
effect  of  the  curvature  in  the  Arrhenius  plot  for  Reaction  (3)  is  quite 
dramatic.  At  low  temperatures,  the  activation  energy  for  this  reaction  is 
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Table  III 


OH+CO— CCU*H 


( 0 K ) 

TEMPERATURE 


( C M3  MOLE  C~3  SEC*) 


kbimolecular* 


2 9 8 

( 1 . 

5 1 1 

( 1 . 

6 0 7 

( 1 . 

7 0 3 

( 2 . 

8 0 0 

( 2 . 

9 2 2 

( 3 . 

9 9 2 

( 3 . 

4 5 ♦ ,1  0 ) x 1 0 _1  5 
7 7 ± .0  8 ) x 1 0 3 
89±.15)xl0'13 

0 3 i .2  6 ) x 1 O'1  3 
39±.12>xl0‘13 
09+.A0)xl0'13 

1 1 1 .4  0 ) x 1 O'1  3 


* Stated  error  bounds  represent  the 


2°  VALUES. 


11 


I 


near  zero;  at  higher  temperatures,  (T ) increases  rapidly  with  temperature. 

[It  is  worth  noting  here  that  visual  comparisons  of  the  'degree'  of  curvature 
in  Arrhenius  plots  for  several  reactions  can  be  misleading;  for  a reaction 
where  the  low- temperature  rate  constant  itself  depends  strongly  on  temperature, 
e.g.,  OH  + CH^  -*•  H^O  + CH^,  substantial  high-temperature  rate  constant  enhance- 
ment may  be  perceived  as  only  moderate  curvature.]  As  in  our  studies  of 
Reactions  (1)  and  (2),  our  measurements  of  ( T)  fill  the  gap  between  500  K 
and  1000  K and  verify  the  adequacy  of  the  empirical  forced  fit  expression 
of  Baulch,  et  al.^  Our  measurements  of  k-j(T)  above  900  K lie  slightly  above 
the  solid  line  in  Figure  5;  it  would  be  worthwhile  to  remeasure  and  extend 
k3(T)  to  the  1200  K temperature  regime. 

Rate  constant  data  for  Reaction  (4),  OH  + C£Hg  -►  H^O  + C^H^.  is  tabulated 
in  Table  TV  and  plotted  in  Figures  6 and  7.  Prior  to  these  measurements,  the 
only  reliable  rate  constants  existing  for  this  reaction  were  a few  measurements 
at  and  slightly  above  room  remperature.  Once  again,  it  can  be  seen  from  Fig- 

g 

ure  7 that  straight-line  extrapolation  of  low-temperature  measurements  re- 
sults in  an  underestimation  of  the  bimolecular  rate  constant  k^(T)  at  higher 
temperatures;  concave  upward  curvature  is  again  evident. 

The  chemical  kinetics  experiments  described  above  represent  valuable  new 
input  to  current  models  of  hydrocarbon  combustion  processes.  Reactions  (1) 
through  (3)  are  of  such  fundamental  importance  to  hydrocarbon  oxidation  chains 
that  accurate  characterization  of  these  k^(T)  functions  is,  in  itself,  a goal 
of  considerable  significance.  Kinetic  estimates  for  Reactions  (1)  - (3)  per- 
vade hydrocarbon  combustion  modeling  efforts;  use  of  incorrect  input  could 
modify  the  bulk  observable  results  of  a simulated  model  process  to  an  extent 
which  implies  that  the  assumed  model  structure  is  wrong  - this  may  be  the  case, 
or  the  problem  could  be  in  the  inaccuracy  of  the  data  base. 


Table  IV 


i 

I 0H+C2H,-H20+C2H5 


I 

(°K) 

( C M ^ HOLEC'^  SEC”-*"  ) 

I 

TEMPERATURE 

kb 

IMOLECULAR* 

I 

T 

297 

(2.59 

± . 21  ) X 

IO-I5 

m 

?• 

400 

(7.71 

± . 76)  X 

1 0 _ 1 3 

i. 

I" 

499 

(1.58 

± . 10)  X 

10‘12 

1* 

»• 

609 

(2.61 

± . 33)  x 

10'12 

1. 

697 

(3.65 

± . 25)  x 

io-12 

li 

800 

(5.07 

± . 34  ) x 

10‘12 

* Stated  error  bounds  represent  the  2°  values. 
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Yet  one  should  not  assume  that  the  value  of  these  results  (and  of  our 
developed  experimental  technique)  is  limited  to  the  numerical  achievements 
detailed  above.  Typically,  the  complexity  of  the  chemistry  involved  in  a 
combustion  and/or  plume  process  is  nearly  overwhelming.  In  recent  years 
it  has  become  increasingly  obvious  that  it  is  simply  not  possible  to  measure 
the  kinetic  parameters  for  all  contributing  elementary  reactions  in  a com- 
plex process.  Indeed,  if  such  a data  base  did  exist,  computational  limita- 
tions would  prevent  its  full  utilization.  Increasing  attention  is  thus 
being  focused  on  the  generation  of  'global'  kinetic  equations,  i.e.,  kinetic 
expressions  which  are  approximately  applicable  to  families  of  analogous 
chemical  reactions.  Such  a viewpoint  then  argues  the  need  for  investigations 
which  emphasize  the  effects  of  chemical  structure  on  the  rates  of  reactivity; 
studies  aimed  at  generating  an  improved  understanding  of  the  general  tempera- 
ture dependence  of  bimolecular  reaction  rates  must  be  pursued.  Investiga- 
tions of  Reactions  ( 1 ) -(4)  and  our  planned  experiments  on  processes  such  as 


OH  + D2  - OHD  + D (5) 

OD  + CO  -*•  D + C02  (6) 

OH  + C3He  - Hz0  + C3H7  (7) 

OH  + CH20  -»■  H20  + CHO  (8) 


r 

i» 


represent  such  structure-motivated  studies. 
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PUBLICATIONS 


"Absolute  Rate  Constant  Measurements  of  the  Reactions  of  OH  with  H2, 

CH4,  and  CO  at  Intermediate  Temperatures,"  F.  P.  Tully  and  A.  R. 

Ravishankara,  J.  Chem.  Phys.  (in  preparation). 

id 

NOTE : Ideally  we  hope  to  group  the  reporting  of  our  data  according  to 
structure,  e.g., 

1)  OH  + H2,  O2  -*•  Products 

2)  OH,  OD  + CO  -*>  Products 

3)  OH  + CH„,  C2Hs,  CsHe  Products 

The  likelihood  of  this  format  depends  on  the  rapidity  with  which  data  for 
Reactions  ( 5 ) - ( 7 ) can  be  obtained. 

PERSONNEL 

Although  most  of  these  experiments  have  been  performed  by  the  Co- 
Principal  Investigators,  undergraduate  student  assistants  have  made 
significant  contributions  to  the  execution  of  this  program.  In  particular, 
students 

K.  Carr  - part-time  assistant  and  Special  Projects 
Researcher  in  Chemistry 

R.  Thompson  - Cooperative  Division  Student 

have  benefited  the  program  and  benefited  from  their  experience. 
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